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Ocean Raman Scattering in Satellite Backscatter UV Measurements 

Alexander P. Vasilkov 1 , Joanna Joiner", James Gleason 2 , Pawan K. Bhartia 2 

Light from the sun, reflected from the Earth’s surface, passing through the atmosphere, is 
measured by remote sensing instruments in space. The light from the sun, the solar spectrum, is 
not smooth. Examined in detail the solar spectrum has a considerable amount of structure, called 
Fraunhofer lines, after Joseph von Fraunhofer who observed these lines in 1817. Trace species 
like ozone and aerosols absorb some of this light and the amount of absorbance is used to infer 
the amounts of these trace species in the atmosphere. The absorbance measurements compare 
direct light from the sun with reflected light from the Earth. 

In the past it was assumed that the spectral structure from the sun, present in both the direct 
and reflected light cancelled out, leaving only the changes induced by the Earth’s surface, 
absorption by gases in the atmosphere, and Rayleigh scattering by air molecules (which produces 
the blue color as seen with human eyes). Previous work from this group showed that this was an 
incorrect assumption and that small changes in the reflected Fraunhofer spectrum were caused by 
a different kind of scattering of the sunlight by molecules in the atmosphere, called Raman 
scattering. Most light undergoes Rayleigh scattering by air molecules. In this kind of scattering, 
the wavelength of light is not changed by the interaction with an air molecule. However, a small 
fraction of light is Raman scattered, which changes the wavelength of the scattered light. Due to 
the Fraunhofer line structure in the incoming solar light, Raman scattering causes the lines in the 
Earth-viewing satellite data to look filled in as compared with the incoming solar light. 

This current work extends that analysis to show that there are additional changes in the 
changes in the reflected Fraunhofer spectrum caused by the Raman scattering of the sunlight by 
water molecules in the ocean. The paper describes the physical basis for the observations, the 
observed spectral changes, and shows how the spectral changes may be quantified to estimate the 
chlorophyll content in the ocean. This approach to measuring ocean chlorophyll complements 
existing methods, because it requires spectrally resolved measurements over a small wavelength 
range and is relatively insensitive to the absolute calibration of a satellite instrument. 
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Abstract 

Ocean Raman scattering significantly contributes to the filling-in of solar Fraunhofer lines 
measured by satellite backscatter ultraviolet (buv) instruments in the cloudless atmosphere over 
clear ocean waters. A model accounting for this effect in buv measurements is developed and 
compared with observations from the Global Ozone Monitoring Experiment (GOME). The 
model extends existing models for ocean Raman scattering to the UV spectral range. Ocean 
Raman scattering radiance is propagated through the atmosphere using a concept of the Lambert 
equivalent reflectivity and an accurate radiative transfer model for Rayleigh scattering. The 
model and observations can be used to evaluate laboratory measurements of pure water 
absorption in the UV. The good agreement between model and observations suggests that buv 
instruments may be useful for estimating chlorophyll content. 

Introduction 

The filling-in and depletion of solar Fraunhofer lines, known as the Ring effect [Grainger and 
Ring, 1962], is a significant component of radiances observed by backscatter ultraviolet (buv) 
satellite instruments. Rotational-Raman scattering (RRS) of atmospheric N 2 and O 2 is a major 
contributor to the Ring effect. Models for RRS have been developed and have compared well 
with observations from the satellite backscatter ultraviolet (SBUV) spectrometer and the Global 
Ozone Monitoring Experiment (GOME) [see e.g. Joiner et al., 1995; Vountas et al., 1998]. 


1 



Properties of RRS can be exploited to retrieve cloud pressure from buv observations [Joiner 
and Bhartia, 1995; de Beek et al. , 2001]. Accurate modeling of the Ring effect is crucial for 
retrieval of cloud pressure as well as trace gas (e.g. NO 2 and SO 2 ) amounts from instruments 
such as GOME and the Ozone Monitoring Instrument (OMI) to fly on NASA’s EOS Aura 
satellite [see e.g. Vountas et al., 1998]. Derived cloud pressures from these instruments may be 
used to improve retrieval of total column ozone and other trace gases. 

Raman scattering in the ocean has been not been included in these studies. Ocean Raman 
scattering has been observed and modeled in the visible spectral range [e.g. Sugihara et al., 
1984; Marshall and Smith, 1990; Haltrin and Kattawar, 1993; Sathyendranath and Platt, 1998] 
and causes filling in of solar Fraunhofer lines. 

Figure 1 shows the ratio of the observed filling-in from SBUV continuous spectral scan 
measurements (taken approximately once per month from 1979-1986) to that computed using the 
Joiner et al. [1995] atmospheric RRS model over ocean scenes with low reflectivity (relatively 
cloud free). A similar spatial pattern with temporal variability is observed in GOME data. The 
filling-in over clear waters can exceed 40% of that expected from the atmosphere. There is 
spatial anti-correlation between the excess filling-in and chlorophyll concentrations derived from 
e.g. (ref). This suggests that the excess filling-in is a result of ocean Raman scattering. 

In the present paper we extend existing models for ocean Raman scattering to the UV spectral 
range and describe a method to account for the effect in satellite buv measurements. The model 
is compared with data from the GOME instrument. The optical properties of water are not well 
known in this spectral range so that comparisons between the model and observation may be 
used to evaluate laboratory measurements of these properties. 
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Figure 1: Excess filling-in at the Ca K line for scenes with reflectivity < 25% from SBUV 
continuous scan data from 1978-1986. 

Ocean Raman Model 

Top-of-atmosphere (TOA) radiance measured by a satellite instrument, /, can be partitioned 
into two components: one consisting of photons that never penetrate into water, I a , and photons 
scattered at least once in water /„., / =/ a +I W , where the subscripts a and w will be used to 

denote atmosphere and ocean, respectively. Following Joiner et al. [1995], we quantify the 
Raman scattering effect on TOA radiance by defining a filling-in factor: 


*(A) = 


/-/' 

/' 


I a +I„-U'a+0 

/'+/! 
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where quoted quantities denote radiances calculated with only elastic scattering. All quantities 
are convolved with an appropriate instrument band pass. The filling-in factor can be expressed as 
a sum of the filling-in factor calculated without Raman scattering in water, k 0 (A), and the 

filling-in excess due to Raman scattering in water, k w (A), i.e. k(A) = k 0 (A) + k w (A) . Here, we 
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compute Jfc 0 (A) using the RRS model of Joiner et al. [1995]. This scheme uses the concept of 
Lambertian equivalent reflectivity (LER) that assumes a Lambertian reflecting surface imitating 
aerosol scattering, Fresnel reflection from the ocean surface, and backscatter in water. We will 
henceforth follow the LER concept in modeling the ocean Raman effect. 

The ocean Raman scattering filling-in factor, k w (X ) , can be expressed as a sum of two terms. 

The first term, k w i(X), is positive and represents energy transferred form shorter excitation 
wavelengths, X e , to an observed wavelength A: 

)_ F„a)p,w f.wp.w ’ ' 

where F 0 is the solar flux, F 0 is the solar flux convolved with the instrument bandpass, and p is 


the TOA reflectance. The integration in (2) is performed over excitation wavelengths of the 
ocean Raman scattering band \Haltrin and Kattawar, 1993]. The second term, k W 2 (X), is negative 
and is due to Raman scattering that transfers energy from the observed wavelength to longer 
wavelengths: k K2 (X) = p*(A',A)/p fl (A). The TOA reflectance, p a (A) , is calculated using the 
LER. The central wavelength of the Raman scattering band is determined by a frequency shift of 
Av=3340 cm 1 . 


Assuming clear skies, the Raman TOA reflectance can be expressed in the form: 


* a A s = nKAm A,e)R„a,A,) 


l-S fc (A e )A 


where T is the atmospheric transmittance, 6b is the solar zenith angle, 6 is the viewing zenith 
angle, R rs is the remote- sensing reflectance just above the ocean surface (Fresnel reflection not 
included), S b is the fraction of solar flux backscattered by the atmosphere toward the ocean 
surface, and A Is the LER. The term 1 - S b (X e )A accounts for an increase of surface irradiance at 
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the excitation wavelength due to multiple reflections between the ocean and the atmosphere. The 
remote sensing reflectance can be expressed through the irradiance reflectance, R, just beneath 
the ocean surface. We assume isotropic angular distribution of the upwelling Raman-scattered 
radiance. 

To calculate the Raman reflectance just beneath the ocean surface, we use the major term of 
the Raman reflectance formula derived by Sathyendranath and Platt [1998]. After some 
transformation this formula becomes 

m, k ) = . < 4 ) 

0 a,) + 6 s a,) + ^^|a(A)+6 > (A)| 

Pu (A) 

where b R is the Raman- backscattering coefficient, a is the absorption coefficient, bb is the 
backscattering coefficient, fid is the mean cosine for the downwelling irradiance, fi R is the mean 
cosine for the upwelling Raman- scattered irradiance. Because of the symmetric form of the 
Raman phase function, /i* ~ 0.5. The quantity /id varies between 1 and 0.5 depending on the 

angular structure of light incident on the ocean surface. We adopt a mean value of 0.75. To 
calculate the Raman reflectance of the ocean we need to specify the inherent optical properties 
(IOP) of seawater: b R , a, and b h . 

The Raman-backscattering coefficient is half of the total Raman scattering coefficient, i.e., 
b * =0.5 b R , owing to the symmetric form of the Raman phase function. The Raman scattering 

coefficient has been found experimentally to vary with the emission wavelength as A' 5 [Sugihara 
et al, 1984]. For excitation in the UV, the spectral dependence of the Raman scattering 
coefficient is b R (A, ) ~ b R ( A s / A 0 )~ 5 , where b R = 2.6e-4 m' 1 , Ao = 488 nm [ Marshall and Smith, 
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1990], The total IOP are the sums of the IOP of pure seawater and scattering and absorbing water 
constituents. 

At present, there is no consensus on the pure water absorption coefficient in the UV. Pope and 
Fry [1997] and Sogandares and Fry [1997] found that the pure water absorption coefficient is 
significantly lower (by a factor of two) than that given by Smith and Baker [1981]). Pure water 
absorbance substantially determines the magnitude of Raman reflectance. A comparison of the 
available measurements of pure water absorption in the UV is shown in Figure 2. Here, we use 
the interpolated data unless specified otherwise. The pure seawater backscattering coefficient is 
taken 
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Figure 2. Available pure water absorption measurements- Q&I: Quickenden and Irvin [1980]; 
P&F: Pope and Fry [1997]; S&F: Sogandares and Fry [1997]; S&B: Smith and Baker [1981]. 

A model of IOPs in the UV is similar to one proposed by Vasilkov et al. [2001]. The model 
was updated by specifying the chlorophyll- specific absorption coefficient as a function of 
chlorophyll concentration. The particulate matter absorption is expressed through chlorophyll 
concentration, C, and the chlorophyll- specific absorption coefficient; a P (X) = Ca* p (C,A) . 
Parameterization of the chlorophyll-specific absorption coefficient in the UV is similar to the one 
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developed in the visible by Bricaud et al. [1995]: a p (C,A) = A(X)C Ba) . The coefficients A(X) 

and B(X) were determined from CalCOFI data sets [G. Mitchell and M. Kahru, private 
communication, 2001]. 

The IOP model contains three input quantities: the dissolved organic matter absorption 
coefficient at a reference wavelength, ao, the particulate matter backscattering coefficient at a 
reference wavelength, bo, and chlorophyll concentration, C. Chlorophyll is a standard product of 
ocean color sensors such as SeaWiFS and MODIS. To determine other quantities, the Case 1 
water model [Morel, 1988] is assumed. According to the model, the DOM absorption at 440 nm 
is 20% of the total absorption of pure seawater and particulate matter pigments. This assumption 
determines the important parameter ao. The backscattering coefficient, which is normally much 
less than the absorption coefficient in the UV, is expressed through the total scattering 
coefficient. A value of the particulate total scattering coefficient at 550 nm was approximated as 
6(550)= 0.3 C 0 ' 62 [Morel, 1988]. Thus, all the input parameters are expressed as a function of a 
single input quantity - chlorophyll concentration. 

Ocean Raman model calculations 

In this section we present some results of calculations of k w as a function of several parameters. 
All the calculations were performed using a triangular band pass with a full-width-half- 
maximum of 0.45 nm which is the approximate resolution of OMI. Spectral dependence of the 
filling in due to both atmospheric, k 0 (X), and oceanic, k w (X), Raman scattering is shown in 
Figure 3. For shorter wavelengths, the ocean filling-in is negative representing a net depletion of 
energy k W 2 > k W }. The opposite is true at longer wavelengths. The spectral signature is similar for 
ko(X ) and kjX). However, the magnitude of k w decreases with X owing to decreasing amounts of 
radiance reaching the surface at excitation wavelengths in the Huggins ozone band. 
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Figure 3: Right: atmospheric Raman filling-in; Left: oceanic Raman filling-in. 


Figure 4 shows that k w decreases with increasing solar zenith angle (0 O ). This dependence is 
different from that of atmospheric RRS that increases with increasing 0 O [Joiner et ai, 1995]. 
The decreasing the oceanic filling-in is explained by reduction of the atmosphere transmittance 
at excitation wavelengths with increasing 0 O . 


OCEAN Co K LINE FILLING IN, PHI=90, REFL=0.1, CHL=0.f 



Figure 4: & K (0 0 ) for 0 = 0, 15, 30, 45, 60, and 70° (corresponding curves from top to bottom). 
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Figure 5 shows that k w diminishes rapidly with increasing chlorophyll. This illustrates that 
there is potential to estimate chlorophyll concentration, especially low amounts, from an 
instrument such as OMI. 

Figure 6 shows that k w significantly depends on the specification of pure water absorbance. In 
the next section, w r e use GOME observations to independently evaluate these laboratory 
measurements of pure water absorption. 



Figure 5 (left): Chlorophyll dependence of k w for 8 0 = 0, 30, 45, 60, 70, 77, 81, 84, 86, and 88° 
(corresponding curves from top to bottom). 

Figure 6 (right): k»{B 0 ) for different estimates of pure water absorbance. Notations as in Figure 3. 

Comparison with GOME observations 

Figure 7 shows model results using climatological chlorophyll concentrations from ? 
compared with spectra from the GOME instrument for clear and turbid waters. The low LER 
reflectivities (12.6% and 6.8%, respectively) indicate relatively cloud-free pixels. In turbid 
waters off the coast of Mexico, model calculations with and without the ocean Raman 
contribution differ little. However, in the clear waters of the southern Pacific, observations show 
a significant difference from model calculations with only atmospheric RRS. Model calculations 
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with zerofclimatological) chlorophyll content slightly overestimate(underestimate) the filling-in. 
The difference between model and observations may be due to a combination of errors in water 
IOPs and/or the fact that the true chlorophyll content may differ from the climatological value 
used here. These calculations were performed using the interpolated water absorbance from Pope 
and Fry [1997]. As Figure 6 indicates, the data are more inconsistent with the measurements of 
Smith and Baker [1981]. 




Figure 7: GOME spectra from orbit 174 on March 24, 1998 and model calculations with and 
without ocean Raman scattering. Left: 28°S, 121°W; Right: 21°N, 108°W. 

Conclusions and Future Work 

Ocean Raman scattering significantly contributes to the total filling-in of solar Fraunhofer lines 
over clear ocean waters. Our developed model agrees well with observations from GOME and 
favors the measurements of pure water absorption by Pope and Fry [1997]. We plan to use the 
model to estimate of chlorophyll concentration from satellite buv observations, This technique 
will utilize the spectral structure of the ocean Raman effect, which is fundiimentally different 
from conventional ocean color algorithms and could be complementary as it has high sensitivity 
at low chlorophyll concentrations. An advantage of this technique is that because it uses high- 
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frequency spectral structure, it is less affected by errors that have a smooth wavelength 
dependence such calibration error. 
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